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Abstract 



Motivated by the recent observations of the peak effect in high- Tc YBCO 

superconductors, we reexamine the origin of this unusual phenomenon. We 

show that the mechanism based on the /j-dependence (nonlocality) of the 

vortex-lattice tilt modulus (744(k) cannot account for the essential feature of 

the peak effect. We propose a scenario in which the peak effect is related 

to the melting of Larkin domains. In our model, the rise of critical current 

with increasing temperature is a result of a crossover from the Larkin pinning 

length to the length scale set by thermally excited free dislocations. 
74.60. Ge, 64.70.Dv 
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About 35 years ago Le Blanc and Little |l| discovered a striking phenomenon, later 
known as the "peak effect" , in a conventional superconductor Nb that the sample can carry 
more supercurrent at a higher temperature (or field) slightly below Bc2{T) where it becomes 
normal. Over the years, the peak effect was found to be ubiquitous in conventional super- 
conductors and it has been observed recently in the high-T^ superconductor YBCO 
crystals |^J^. Pippard and Larkin and Ovchinnikov p| proposed that the peak effect is 
a result of an anomalous softening of the vortex lattice. The basic physics is that a soft 
lattice can be pinned more strongly than a more rigid lattice. In fact, an infinitely rigid 
lattice cannot be pinned at all by random pinning. The unresolved problem, however, is the 
mechanism which leads to the abrupt loss of the vortex-lattice rigidity. 

In this paper, we first reexamine the standard interpretations of the peak effect and show 
that the mechanism of an anomalous softening of the wavevector dependent tilt modulus 
C44(k) does not account for the essential features of the peak effect. We then propose a 
scenario in which a melting of the "Larkin domains" leads to the peak effect. The rise of the 
critical current with increasing temperature is a result of a crossover of two length scales: 
from the Larkin pinning length to the average separation between thermally excited free 
dislocations 0. 

Let us first recall briefly the general features of the peak effect phenomenon in both con- 
ventional superconductors and high-T^ superconducting YBCO single crystals. In samples 
with high values of the critical current density jc at a fixed magnetic field B, one usually 
finds that jc decreases to zero monotonically with increasing temperature and the peak effect 
is absent But in samples with low jc (weak pinning), the temperature dependence 

of jc can be quite different. Fig. |l] is a plot of the critical current density as a function of 
temperature extracted from Ref. for a YBCO crystal. With increasing temperature, jc 
initially decreases monotonically, then suddenly rises, reaches a peak before finally dropping 
to zero. Experimentally, the peak effect is identified as a dip in resistance [Q-^J^, a dip in 
the in-phase part of ac susceptibility or a peak in critical current density jc (obtained 
with the standard voltage criterion) as a function of T or 5. In Nb and other low-Tc 
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superconductors, the onset (the rise of jc) temperature Tp of the peak effect is very close 
to T^{B) and T^{B) - Tp ^ 0.5 K. In YBCO, T^{B) - Tp ~ 5 K is about 10 times larger. 
However, {Tc{B) — Tp)/Tc{B) ~ 0.05-0.1 is about the same for both low- Tc and high-Tc 
superconductors. The problem under consideration here is why the critical current rises 
with increasing temperature. 

Pippard proposed that the rise of the critical current has to result from a rapidly 
decreasing rigidity of the vortex lattice. The rigidity of the vortex lattice prevents the 
vortex lines from taking advantage of the valleys of random pinning potential. Thus a 
rapidly decreasing rigidity would allow the lattice to conform better to the pinning potential 
and enhance critical current. This idea on pinning was subsequently put forward more 
rigorously by Larkin and Ovchinnikov (LO) in their theory of collective pinning P). It was 



shown by Larkin |T0| that in the presence of random pinning the vortex lattice loses its 
long-range translational order and breaks up into domains of correlated regions in which the 
vortex lines interact elastically. LO argued p] that the critical current density is determined 
by the fluctuations of random potential in a domain and the pinning force density jcB = 
{nf^/VcY^^, "where n and / are the density and strength of the pins, VJ, the volume of the 
domain. The size of the Larkin domains can be estimated by a simple energy consideration. 
The vortex lattice deforms to take advantage of the random pinning potential at the cost of 
the elastic energy. The total unit volume energy change is |§] 

SF = C,,qr + C,,qr-frp{^Y/\ (1) 

where Cqq is the shear modulus of the lattice, Vp the range of the pinning potential, R 
and L are the transverse (to the field) and longitudinal (along the field) dimensions of the 
domain, and V = B?L. The minimization of Eq. (|1|) gives the pinning lengths Rc and L^: 
Rc ~ CgQ^Cll'^rp/np , Lc = {C^^^/C^qY^'^Rc- In very thin samples with a perpendicular 
field, if the pinning is so weak that Lc is greater than the sample thickness, the problem 
becomes two-dimensional (2D) and only Rc ~ C^'^plrt}!'^ f is relevant. In the LO theory, 
since jcB = {nf'^/VcY^'^, ^^e peak effect can be accounted for if the volume of Larkin domain 



Vc drops faster than n/^ in some field or temperature range. The central question here is 
what mechanism does that. 



It was found by Brandt |TT| that, near upper critical field Bc2 the vortex fields overlap 
strongly and the tilt modulus C44 becomes nonlocal: it softens substantially for short wave- 
length tilt deformation. Most experiments of the peak effect on low-T^ materials are carried 
out as a function of field while keeping the temperature constant, the peak effect manifests 
itself as a peak in jc (or a dip in resistance) near Bc2{T). It is thus natural to relate the 
C44(k) softening mechanism to the peak effect. LO found p| that this mechanism leads to 
an exponential form for Vc when Rc becomes smaller than A' = A/(l — bY^'^, where A is 
the penetration depth and b = B/Bc2, and jc ~ exp(— _BCgg^/c/irp/l^), with W = np and 
kh = 1/A'. It was customary to assume a scaling function for the field dependence of W, 
W b{l — by, and with Cqq ~ 6(1 — 6)^ and A;^ ~ (1 — 6), one indeed finds that jc rises 
exponentially with B for 6 ~ 1. 

The above interpretation of the peak effect has two major difficulties. The first is that if 
this mechanism is the relevant one, it should also account for the temperature dependence: 
jc rises with increasing T in the peak regime. Giving the most liberal estimate for the 
temperature dependence of the parameters, however, the above mechanism fails to explain 
why jc rises with T. The elementary pinning force / is a function of the local gradient of 
the amplitude of the gap function, or / ~ ai|A(r)p. For Tc smearing pins, ai ~ (1 — t), 
|A(r)|2 ~ (1 - t), t = T/Tc, / ~ (1 - ty, and for pins that do not smear T^, / ~ (1 - t). 
Thus W {l-tf 01 {l-t)\ Without mehing, ~ (1 - t) and 4 ~ (1 - ty/^. Thus jc 
either does not change with T or decreases exponentially with increasing temperature and 
jc never increases with increasing T, according to this mechanism. The second difficulty of 
this mechanism is that the peak effect has been observed in thin films and in very thin 
NbSe2 crystals with pinning weak enough such that Lc exceeds the sample thickness P,pr3 



in which C44 does not seem to play any role. Therefore we believe that this mechanism 
cannot be responsible for the peak effect. 

Another modulus of the vortex lattice which enters the pinning problem is the shear 



modulus Cee- The peak effect, the rise of jc, suggests to us that Cqq of the vortex lattice 
vanishes at the peak effect regime much faster than what was calculated [|12| from the 
Ginzburg-Landau theory. In fact, Pippard ||^ suggested that the peak effect could be due 
to the softening of the shear modulus of the vortex lattice. Here we propose that the rise of 
the critical current, or the peak effect itself, is a result of vortex lattice melting. 

Indeed, the fact that a vortex lattice can melt has been pointed out by several authors 



1^-|16|. Most of the efforts has been focused on the possibility of melting of a perfect 
Abrikosov lattice and the determination of the phase diagrams. A 2D vortex lattice can 
melt via the mechanism of the unbinding of the dislocation pairs n,|l5|, similar to the 



melting of a 2D crystal studied by Kosterlitz and Thouless [17 , and by Halperin and Nelson 



T8| . Much less is known for the melting of a 3D lattice in general. For a perfect 3D vortex 



lattice both analytic considerations [|T9| and numerical simulations [|20| suggest a first order 
transition. In particular, it has been shown that a finite density of free edge dislocations 
would result in a zero long wavelength shear modulus [ pT| . 

Generally speaking, the melting of a lattice (2D and 3D) in the presence of quenched 
disorder is not well understood. Indeed, often it is not clear what "melting" means in this 
situation. Naively, if the pins are dilute compared to the vortex density, the melting would be 
more or less the same as in a pure system and would result in a free (unpinned) vortex fluid. 
It is doubtful that such a regime exists in real samples. In fact, almost all weak-pinning (low 
jc) samples show peak effect which implies a high density of pins according to our 

scenario described below. It was argued that a vortex lattice, disordered due to random 



pinning, could be a "vortex glass" at low temperatures, and the melting transition is then 
a vortex-glass-to-vortex-fiuid transition, which is second order at least for strong random 
pinning. While a vortex glass phase, if it exists, would determine the small current behavior 
especially for systems with strong disorder, the peak effect is a property of (relatively) large 
current in samples with weak pinning. 

For a vortex lattice in the presence of random pinning, the lattice deforms locally to take 
advantage of the fluctuations of the random pins and the long range order of the lattice 



is destroyed |T^. The lattice is ordered within Larkin domains whose size is the pinning 
lengths Rc and Lc- (We assume weak pinning by which we mean that the typical pinning 
force / for individual pins is small. We shall always assume that the density of pins n is not 
small compared with the vortex density, which seems to be the case even in samples with 
extremely low jc-) Note that Rc and are also the length scales beyond which the lattice is 
elastically decoupled. In other words, a small local shear deformation of the lattice would not 
propagate (elastically) much farther than Rc. How would such a weakly pinned lattice melt? 
Melting of a lattice is usually characterized by the vanishing of the long wavelength shear 
modulus. For a pinned vortex lattice, however, a long wavelength shear is not sustained 
by the long wavelength elasticity of the lattice. If we shear slightly the opposite sides of a 
pinned vortex lattice, the shear deformation would decay (presumably exponentially) inside 
the lattice with a decay length of order Rc, and the vortices inside would not feel the shear. 
In other words, the long wavelength shear is sustained by the pinning force. It is then not 
surprising that the length scale Rc should play a crucial role in the melting of a pinned 
lattice. 

We first consider the case of a 2D lattice where much is known for the melting of a 
clean (pinning free) system |]I^|TB[. In the case of a 2D lattice, thermally excited dislocation 



pairs are bound for temperatures below the melting temperature and, consequently, the 
shear modulus is finite. At T^, the largest dislocation pairs start to dissociate and the 
long wavelength shear modulus drops discontinuously to zero. Above T^, the density of 
free dislocations rises from zero, and the mean distance between free dislocations is the 
Kosterlitz-Thouless correlation length |2^ ( ~ exp[c/(T — T^Y] with u ~ 0.37 |jl8[, which 



diverges as T ^ T+. The correlation length ( also sets the length scale for the g(wavevector)- 
dependence of the shear modulus fi{q,T): roughly speaking, fi is zero for q <1/C and finite 
for q > Now imagine that the vortex lattice is weakly pinned. For T < T^, the Larkin 
length Rc sets the elastic length scale and the critical current density jc = n^^"^ f / RcB. 
At T = T^, the vortex lattice melts with the long wavelength shear modulus dropping to 
zero. However, as argued in the previous paragraph, the pinned lattice would not feel being 



melted at this point since the lattice is elastically decoupled beyond the length scale of Re- 
The critical current density is still determined by Rc- For T > T^, another length scale 
C, the Kosterlitz-Thouless correlation length, enters the system. C, decreases exponentially 
fast from the infinity as the temperature is increased and will soon become comparable 
to Rc- For C < Rci the Larkin domains melt and the relevant elastic length scale for the 
determination of jc is now C,: 

1/2 f 1/2 f 

Thus the onset the of peak effect occurs when the two length scales Rc and ( cross each other 
(Fig. H). The exponential increase of jc with T (Eq. (§)) would continue until ( is of the 
order of the lattice constant. However, the temperature dependence of the pinning force / 
and thermally activated vortex motion [0] would presumably dominate the behavior of jc 
at even higher T and cause jc to vanish. 

In the 3D case, melting of a pure lattice is much less understood. If the melting transition 
is mediated by generation of the free edge dislocations, one would expect a similar mechanism 
for the peak effect as in 2D, with ( now being the mean distance between dislocation lines. 
To make a qualitative or semi-quantitative estimate for jc in the peak effect regime, we take 



the Landau- Ginzburg-like free energy often used in 3D dislocation systems [55 



F{p) = -F,p In Cp + F2P + F^p\ (3) 

where p is the areal density of dislocation lines. Fx and F3 are positive constants, C is a 
constant of the order with a being the lattice constant, F2 > at low temperatures and 
-F2 < at high temperatures. It is easy to see that Eq. (§) implies a first order transition: 

P=\ (4) 

[p,exp[A(T-Tj-], (T>T^) 

where pc = F1/2F3 pa/327T^kBTm, with p being the shear modulus, if we use the values 



in Ref. |]2S| for Fi and F3), A = y— Fg/Fapc, and u = 1/2. The mean distance between 
dislocation lines, ( = p~^/^ = p~^/^ exp[— A(T — TmY is the length scale to be compared 



with the Larkin length R^. In the region where ( < Rc, the critical current density is 
determined by (: 

„l/2 f f 4^ ^ ^.f/ 3 , (( > 

' ^'■'=s"liP^«^e.p|i^,r-r„)1,(c<«.) 

where in the region of C < -Rc Eq- (ID is minimized, with R^ replaced by C, to determine L^. 
Since C, has a discontinuous jump at (if the transition is first order), it is possible that 
jc will have a jump at the onset of the peak effect which may occur in samples with "very 
weak" pinning (Fig. 0). In fact, a jump in was observed experimentally [^. Note that 
our argument does not depend on the detailed nature of the melting transition, e.g. first 
vs. second order |^^, although we have used Eq. (^ to obtain some estimates. The key 



point in our scenario is that some other (elastically relevant) length scale enters the system 
around the melting transition and it crosses the Larkin length. In 2D, we believe this length 
scale is the Kosterlitz-Thouless correlation length. While in 3D, it is most likely that it is 
the length scale set by thermally excited edge dislocations. 

In conclusion, a quantitative analysis suggests that the non-local effect of C44 is not the 
cause for to rise with increasing temperature at the peak effect regime. We propose a 
scenario based on the melting of the vortex lattice in the presence of weak pinning. The 
onset of the peak effect is the crossover of the two elastically relevant length scales. We 
believe this captures the basic physical picture, at least for weak enough pinning, although 
ideally the analysis should be put on a more rigorous footing. On the other hand, if the 
pinning is so strong that the Larkin length is of the order of the lattice constant, the sample 
should not, according to our scenario, show the peak effect. In this paper we have focused 
only on the behavior of jc- The dynamics above jc can be quite different in the peak regime, 
where a significant enhancement of plastic flow was recently observed in 2H-NbSe2. This 
indicates that the dynamically generated defects proliferate in the peak regime as well. We 
emphasize the difference between the scenario concerning melting described in this paper 
and that suggested in literature in which the onset of dissipation in the presence of a 
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driving force or the vanishing of jc has been interpreted as the evidence for the vortex-lattice 
melting transition. The rationale behind this suggestion is based upon the assumption that 
when the vortex lattice melts, the vortex lines will flow. This scenario may be correct if the 
sample is so pure that the vortex line density is much greater than the impurity density. 
Unfortunately, this cannot be achieved even in the cleanest crystals of 2H-NbSe2 where the 
critical current density is many orders of magnitude smaller than that of YBCO. Instead, 
the clean, high quality crystals of 2H-NbSe2 always exhibit |Q a pronounced peak effect and 
high-quality YBCO crystals are also found recently |^,|[| to show a peak effect. 

Note added: In the course of writing this manuscript, we received a preprint by Larkin, 



Marchetti, and Vinokur |3^, in which they attribute the peak effect to the softening of Cqq 
just below the melting transition. 
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FIGURES 

FIG. 1. The temperature dependence of critical current density for a YBCO crystal in a 
magnetic field, extracted from Ref. [5]. We define Tp as the onset temperature at which starts 
to increase. 

FIG. 2. Schematic behavior of the two length scales Rc and C,, as functions of temperature: 
(a) two dimensions; (b) three dimensions with weak pinning; and (c) three dimensions with "very 
weak" pinning. The temperature for the onset of the peak effect is Tp, where C becomes smaller 
than Rc- In (c), Tp = T^. The solid part of the lines determines the critical current density. 
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Figure 1 (Ling, Tang, Bhattacharya, Chaikin) 




Figure 2 

(Ling, Tang, Bhattacharya, Chaikin) 



